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ABSTRACT 

We probe the role of carbon in the nltraviolet (UV) extinction by examining 
the relations between the amonnt of carbon reqnired to be locked np in dnst 
[C/H]j^g^ with the 2175 A extinction bnmp and the far-UV extinction rise, based 
on an analysis of the extinction cnrves along 16 Galactic sightlines for which the 
gas-phase carbon abundance is known and the 2175 A extinction bump exhibits 
variable strengths and widths. We derive from the Kramers-Kronig 

relation which relates the wavelength-integrated extinction to the total dust vol¬ 
ume. This approach is less model-dependent since it does not require the knowl¬ 
edge of the detailed optical properties and size distribution of the dust. We also 
derive [C/H]^^g^ from htting the observed UV/optical/near-infrared extinction 
with a mixture of amorphous silicate and graphite. We hnd that the carbon 
depletion [C/H]^^g^ tends to correlate with the strength of the 2175 A bump, 
while the abundance of silicon depleted in dust shows no correlation with the 
2175 A bump. This supports graphite or polycyclic aromatic hydrocarbon (PAH) 
molecules as the possible carrier of the 2175 A bump. We also see that 
shows a trend of correlating with l/i?v, where Ry is the total-to-selective ex¬ 
tinction ratio, suggesting that the far-UV extinction is more likely produced by 
small carbon dust than by small silicate dust. 

Subject headings: dust, extinction — ISM: abundances — ISM: clouds 


1. Introduction 

Dust is a ubiquitous feature of the cosmos. It has dramatic effects on the interstel¬ 
lar medium (ISM) by affecting the interstellar conditions, participating in the physical and 
chemical processes occuring in the ISM, and shaping the appearance of dusty objects as an 
absorber, scatterer, and emitter of electromagnetic radiation, and as an efficient catalyst for 
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the formation of H 2 and an agent for photoelectrically heating the interstellar gas. To under¬ 
stand the interaction of the solid dust particles in the space between stars (i.e., interstellar 
dust) with radiation and the dust-related interstellar processes, the knowledge of the dust 
composition and size is required. 

What is interstellar dust made of? The ubiquitous detection of the 9.7 pm Si-0 stretch¬ 
ing and 18 pm O-Si-0 bending absorption features in interstellar regions implies that silicate 
must be a major component (Henning 2010). These features are broad, smooth and feature¬ 
less, suggesting that interstellar silicate is amorphous (Li & Draine 2001a, Kemper et al. 
2004, Li et al. 2008). Another absorption feature ubiquitously seen in the diffuse ISM of 
the Milky Way and external galaxies is the 3.4 pm feature, commonly attributed to the C- 
H stretch in some sorts of aliphatic hydrocarbon materials (see Pendleton & Allamandola 
2002). Depending on their H fractions and the preparation method, the experimentally- 
synthesized aliphatic hydrocarbon analogues show considerable variation in the C-H band 
strength (e.g., see Furton, Laiho, & Witt 1999). Therefore, it has not been possible to use 
the measured 3.4 pm feature to determine the abundance of aliphatic hydrocarbon material 
in the ISM (e.g., see Chiar et al. 2013). 

One may infer the amount of silicate dust Mgii from the interstellar Si, Mg, and Fe 
abundances. Let [X/H]jgj^ be the total interstellar abundance of element X relative to H, 
[X/H]g^g be the amount of X in the gas phase, and [X/Hj^^g^. be the amount of X contained in 
dust. As an element will be in gas-phase and/or locked up in solid-phase, we will obviously 
have [X/H]^^g^ = [X/Hjjgjyj — [X/Hj^g^g. Let px be the atomic weight of X (px = 12, 16, 56, 
24, and 28 for carbon, oxygen, iron, magnesium, and silicon). If we assume a stoichiomet¬ 
ric composition of Mg 2 xFe 2 (i_x)Si 04 for silicate (i.e., each silicon atom corresponds to four 
oxygen atoms), the interstellar silicate dust mass relative to hydrogen is 

Msu/Mh = pFe [Fe/H]j^^, + pMg [Mg/H]^^^^ + pa [Si/H]^^^^ + 4 x po [Si/H]^^^^ . (1) 

The exact interstellar abundances for the dust-forming elements are unknown (see Snow & 
Witt 1996, Sofia 2004, Li 2005, Jenkins 2009). The abundances of the Sun, the proto-Sun, 
unevolved early B stars, and F and G stars are often adopted as the interstellar reference 
abundance standards. If we adopt the proto-Sun abundances of [Fe/Hj^ = 34.7 ± 2.4 ppm, 
[Mg/Hj^ = 41.7 ± 1.9 ppm, and [Si/H]Q = 40.7 ± 1.9 ppm (Fodders 2003) for the interstellar 
reference abundances and assume that Fe, Mg and Si are essentially completely locked up 
in dust (i.e., [X/Hj^^^^ = [X/H]jgj^ for Fe, Mg and Si), we obtain Msh/Mh ~ 6.69 x 10“^. 
Let Gabs(9.7pm)/I/ = 1.00 x lO^cm”^ be the peak absorption cross section per volume of 
silicate dust at 9.7 pm (Draine & Lee 1984) and Ay/N-n = 5.3 x 10“^^ magcm^ be the 
visual extinction Ay per hydrogen column A^h (Whittet 2003). One may estimate the ratio 
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of the optical depth of the 9.7 /xm silicate absorption feature to Ay to be 

ATg,T/Av ^ [Cabs(9.7/im)/l/] /psil X (Msii/Mn) x (Nn/Av) x mn ^ 1/16.7 , (2) 

where psii ~ 3.5 g cm“^ is the mass density of silicate and mn = 1.66 x 10“^^ g is the mass of a 
hydrogen atom. This is close to the mean value of Argj/Av ~ 1/18.5 observed for the local 
diffuse ISM (see Draine 2003a), suggesting that the proto-Sun Fe, Mg and Si abundances 
are sufficient to account for the observed 9.7 pm silicate absorption featureji] 

However, silicate alone is not sufficient to account for the observed interstellar extinction. 
As shown by Purcell (1969), the Kramers-Kronig relation can be used to constrain the dust 
quantity in the ISM. Let Ax be the extinction at wavelength A and Ax/Nu be the extinction 
per H nucleon. As shown in Li (2005), the Kramers-Kronig relation relates the wavelength- 
integrated extinction Ax/Ni^dX) to the total volume occupied by dust per H nucleon 

(Kdust/H) through 

[ ^dX = 1.086 X 37r^F%^ , (3) 

Jo Ah H 

where F is a dimensionless factor and depends only upon the grain shape and the static 
(zero-frequency) dielectric constant Eq of the grain material (see below). The interstellar ex¬ 
tinction per H nucleon (Aa/Ah) is not known for all wavelengths, but for a limited range of 
wavelengths. We approximate J^Ax/NudXhy Aa/Ah dA, with the latter derived 

from the silicate-graphite-PAH model of Weingartner & Draine (2001; hereafter WDOl) 
which closely reproduces the observed extinction from the far-UV to the near-IR and up to 
A ~ 30pm (see Figure 16 of Li & Draine 2001b). The WDOl model gives a reasonably ac¬ 
curate estimate for the extinction at A > 30 pm which is observationally unknown since this 
model closely hts the observed dust thermal emission from the near-IR to submillimeterj^ 


^ If we adopt the abundances of unevolved early B-type stars ([Fe/H]^ = 27.5 ± 2.5 ppm, [Mg/H]^ = 
36.3 ±4.2 ppm, and [Si/H]^ = 31.6 ± 1.5 ppm; Przybilla et al. 2008) for the interstellar reference abundances, 
we obtain Atq.t/Av ss 1/21.0. This is somewhat lower than the observed value of Arg.y/Ay Ri 1/18.5. The 
abundances of the silicate-forming elements of the Sun ([Fe/Hjg = 31.6±2.9 ppm, [Mg/H]g = 39.8±3.7 ppm, 
and [Si/H]g = 32.4 ± 2.2 ppm; Asplund et al. 2009) and of young F and G stars ([Fe/H]^ = 28.2 ± 7.8 ppm, 
[Mg/H]^ = 42.7 ± 16.7 ppm, and [Si/H]^ = 39.8 ± 12.8 ppm; Sofia & Meyer 2001) are somewhat intermediate 
between that of the proto-Sun and early B stars, except the magnesium abundance of young F and G stars 
is a little bit higher than that of the proto-Sun and early B stars (but we note that the abundances of young 
F and G stars have large uncertainties; see Sofia & Meyer 2001). 

^The WDOl model reproduces the observed UV/optical extinction as well as the albedo (see Figure 15 
of Li & Draine 2001b). This implies that the WDOl model gives a reasonably accurate description on how 
interstellar dust absorbs starlight. On the other hand, the fact that the WDOl model closely reproduces 
the observed IR emission indicates that this model also gives a reasonably accurate description on how 
interstellar dust emits in the IR and of its IR opacities. 
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For the diffuse ISM, the WDOl model gives Aa/A^h dX ~ 1.49 x 10 magcm^ H 

■^912A 

Since Ax is a positive number and Vdust/A^n is obtained by integrating Ax/N^ over all wave¬ 
lengths (see eq.E]), therefore, the integration of Ax/N^i over a finite wavelength range can be 
used to place a lower limit on the total amount of dust in the ISM. 

If we take the dust to be spheroids with semiaxes a, 6, 6 (prolate if a/h > 1, oblate 
if a/h < 1), then F is related to the static dielectric constant eo cind the “depolarization 
factors” La and L;, = (1 — La)/2 through (Draine 2003b) 


F(a/6,£o) 

where for prolates 

and for oblates 

where 


(^0 - 1 ) 
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e = jl - . 


(4) 

(5) 

( 6 ) 
(7) 


In Figured] we show the F factors for silicate, Fe 203 , amorphous carbon, and graphite. 
These dust species are often invoked to account for the interstellar extinction. They also 
represent materials from being dielectric (silicate with Sq ~ 6.40 Fe 203 with Sq ~ 16) to 
somewhat conducting (amorphous carbon with eg ~ 160) to highly conducting (graphite 
with Sq —)■ cxo). For silicate spheres we have F ^ 0.64. For modestly elongated or flattened 
silicate spheroids, F does not exceed unity: F ^ 0.73 for a/b = 3 (Greenberg & Li 1996)@ 


^In the far-IR, the extinction rapidly declines with wavelength (e.g., Ax oc A“^, see Table 6 of Li & 
Draine 2001b) and therefore, the contribution of the extinction at A > 1000 to the extinction integration 
is negligible. Indeed, the wavelength-integrated extinction f Ax/N^dX only increases ^^1.6% by extending 
the upper integral limit from A = 1000 /rm to 10,000 p,m. 

'^Davis et al. (1988) measured the temperature dependence of the static dielectric constant Eq of natu¬ 
rally occurring monocrystalline forsterite Mgi. 8 eFeo.i 4 Si 04 . They found that Sq depends linearly upon the 
temperature: Eq ~ 6.3 -I- 0.0064 x T. For typical dust temperatures of T ss 20K in the ISM (Li & Draine 
2001b), one obtains Eq ~ 6.4 for silicates. 

^Greenberg & Li (1996) modeled the 9.7 and 18 p,m silicate polarization features toward the Becklin- 
Neugebauer (BN) object in terms of core-mantle spheroids. They found that the a/b = 3 prolate shape 
provides the best fit. 
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and F ^ 0.68 for a/b =1/2 (Lee & Draine 1985)0 Assuming the proto-Sun abundances for 
the interstellar reference abundances, the total silicate volume per H nucleon is 

l/sii/H (Msii/Mn) /psii 3.17 x 10"^^ cm^ . (8) 

With F 0.7 for moderately elongated or flattened silicate spheroids, the wavelength- 
integrated extinction arising from the silicate dust component is at most ~ 7.14x 10“^® mag cm^ 
Apparently, this is not sufficient to account for the observed lower limit of ~ 1.49 x 10“^® mag cm^ H“^. 
Therefore, there must exist other dust component (s), with carbon dust being the most pop¬ 
ular onejfl As a matter of fact, all modern dust models invoke some form of carbonaceous 
materials including graphite (Mathis et ah 1977, hereafter MRN; Draine & Lee 1984), “or¬ 
ganic refractory” (Li & Greenberg 1997), and hydrogenated amorphous carbon (HAG; Duley 
et ah 1989, Jones et ah 2013). We note that carbon-bearing dust species have revealed their 
presence in the ISM through the 3.4/im aliphatic G-H absorption feature (Pendleton & 
Allamandola 2002), the “unidentihed” infrared emission (UIE) bands at 3.3, 6.2, 7.7, 8.6, 

11.3 and 12.7/im (Leger & Puget 1984, Allamandola et al. 1985), and the “extended Red 
emission” (ERE) attributed to the photoluminescence of some form of carbon dust (Witt & 

Vijh 2004). 

What role does the carbonaceous dust component play in the UV/optical extinction? as 
elaborated above, it must account for ~ 52% of the total interstellar extinction. Spectroscop¬ 
ically, the 2175 A extinction bump, the strongest absorption feature observed in the ISM, 
was widely attributed to small graphite grains (Stecher Sz Bonn 1965, Draine & Malhotra 
1993, Mathis 1994) or polycyclic aromatic hydrocarbon (PAH) molecules (Joblin et al. 1992, 

Li & Draine 2001b, Gecchi-Pestellini et al. 2008, Steglich et al. 2010, Mulas et al. 2013). 

It was thought that the tt-tt* transition in hybridization of aromatic carbon causes the 


®The a/b = 1/2 oblate shape was shown to provide a better match than any other shapes to the 3.1 ice 
polarization feature (Lee & Draine 1985) and the BN 9.7 /rm silicate polarization (Hildebrand & Dragovan 
1995). Unlike Greenberg & Li (1996), Hildebrand & Dragovan (1995) considered bare silicate dust and only 
modeled the 9.7 /rm polarization. 

^This reduces to w 5.68 x 10“^®magcm^ H“^ if the Fe, Mg and Si abundances of unevolved early B stars 
(Przybilla et al. 2008) are adopted. The Rn/H values obtained with the abundances of young F and G stars 
(Sofia & Meyer 2001) or the Sun (Asplund et al. 2009) are intermediate between that of the proto-Sun and 
that of the early B stars. 

®Other dust species such as iron and iron oxides may also be present in the ISM. However, if a substantially 
large fraction of the interstellar Fe/H abundance is consumed by iron or iron oxides, the interstellar silicate 
grains will not be “dirty” enough to account for the observed optical/near-IR extinction (see Draine & Lee 
1984). Moreover, iron oxides have absorption peaks at ^9.2, 18, 21 and 30/rm for hematite a-Fe 203 and 
^ 17 and 25/xm for magnetite Fe 304 which are not seen in the ISM (see Zhang, Jiang, & Li 2009). 
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2175 A extinction bump. However, Parvathi et al. (2012) recently examined the correlation 
of the carbon depletion (i.e., the number of carbon atoms [relative to H] locked up 

in dust) with the strength of the 2175 A extinction bump for 21 Galactic lines of sight and 
found no correlation. This provides a severe challenge to the hypothesis of graphite or PAHs 
as the carriers of the 2175 A bump. 

Parvathi et al. (2012) derived the dust-phase carbon abundance for each sight- 

line by subtracting the gas-phase carbon abundance from the interstellar reference 

abundance [C/H]jgj^: = [C/H]jgj^ — [C/H]gg^g. The gas-phase carbon abundance 

[C/H]gg^ was determined by Soha et al. (2011) for six translucent sightlines and by Parvathi 
et al. (2012) for 15 sightlines. Their analysis was based on the strong 1334 A absorption line 
of CII observed by the Space Telescope Imaging Spectrograph (STIS) on board the Hubble 
Space Telescope (HST). As described in Parvathi et al. (2012), these sightlines span a variety 
of Galactic disk environments which are characterized by different extinction and sample 
paths ranging over three orders of magnitude in (rin), the average density of hydrogen. 
They have all been observed with the Far-Ultraviolet Spectroscopic Explorer (FUSE) and 
had their total neutral hydrogen (i.e., HI and H 2 ) column densities measured by Gartledge 
et al. (2004, 2006). Parvathi et al. (2012) took the total interstellar carbon abundance 
to be [G/H]jgj^ = 464 ppm. They assumed that the gas-phase carbon abundance toward 
HD 206773 ([G/Hj^g^g 464 ± 57ppm) represents the lower limit for the interstellar cosmic 
carbon abundance. 

In this work we revisit the Galactic sightlines considered by Parvathi et al. (2012) 
to probe the role of the carbonaceous dust component in the UV extinction, with special 
attention paid to the 2175 A extinction bump. We take two approaches both of which differ 
from that of Parvathi et al. (2012). One approach is independent of any dust models: we use 
the Kramers-Kronig relation of Purcell (1969) to relate the wavelength-integrated extinction 
to the total dust volume (see eq.E]) to infer the carbon depletion [G/Hj^^^^ for each sightline 
(see [|2]). We then examine the relation between [G/H]^^^^ and the 2175 A extinction bump 
strength and the far-UV extinction rise (see [|2]). Another approach is to model the observed 
extinction curve of each sightline in terms of the silicate-graphite model (see [|3]). Again, 
the required carbon depletion in graphite is compared with the 2175 A extinction bump 
strength and the far-UV extinction rise (see [12]) • The results are discussed in [g] and the 
major conclusions are summarized in [|5l 
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2. Carbon Depletion Inferred from the Kramers-Kronig Relation 

Among the 15 Galactic sightlines studied by Parvathi et ah (2012) and the six sight- 
lines studied by Soha et ah (2011), we select 16 sightlines for which the extinction curves, 
hydrogen column densities, and gas-phase carbon abundances have been well determined. 
In Figures mE] we show the UV/optical/near-IR extinction curves of these 16 sources con¬ 
structed as following, with the extinction parameters taken from Valencic et ah (2004) and 
Gordon et ah (2009). 

The UV extinction at 3.3 < < 8.7/rm”^ is derived from the following analytical 

formula which consists of three parts: a linear background, a Drude profile for the 2175 A 
extinction bump, and a far-UV nonlinear rise at > 5.9 /rm“^ (see Valencic et al. 2004, 
Gordon et al. 2009): 


Ax!Ay = c'l -f 411(0:, 7 ,Xo) + c!iF{x) 


(9) 


F(x) 


Zl(x,7,Xo) 


(x 2 


x^ 

Xq)^ -1- X272 


0 , 

0.5392 (x - 5.9)^ + 0.05644 (x - 5.9)^ 


X < 5.9 /im 
X > 5.9 /im 


( 10 ) 

( 11 ) 


where x = 1/A is the inverse wavelength in /im“^, 4 and 4 define the linear background, 
4 defines the strength of the 2175 A extinction bump which is approximated by Zl(x,7,Xo), 
a Drude function which peaks at Xq ~ 4.6 /im“^ and has a FWHM of 7 , and 4 defines 
the nonlinear far-UV rise. This parametrization was originally introduced by Fitzpatrick & 
Massa (1990; hereafter FM90) for the interstellar reddening 


E{X-V)/E{B-V) = Rv{Ax/Av-l) = ci + C2X + c^D{x,'^,Xo) + CiF{x) , ( 12 ) 


where E{\ — V) = Ax — Ay, E{B — V)= As — Ay, and As is the R-band extinction. It is 
easy to hnd the relationship between the FM90 parameters and that of Valencic et al. (2004) 
and Gordon et al. (2009): 


,^(cj/Rv + l, J = l , 

\c,/Ry, j = 2,3,4 , ^ ^ 

where Ry = Ay /E{B — V) is the total-to-selective extinction ratio. The c' parameters of the 
selected 16 sources are taken from Valencic et al. (2004) and Gordon et al. (2009). Finally, 
for a given Ry, the optical/near-IR extinction at 0.3 < A“^ < 3.3 /im“^ is computed from the 
GGM parametrization (see Gardelli et al. 1989). We then smoothly join the UV extinction 
at A“^ > 3.3/im“^ to the optical/near-IR extinction. 



Let Aint = Aa/A^h (iA be the extinction integrated over all wavelengths. The total 
volume of carbon dust per H nucleon (Vc/H) can be derived from the Kramers-Kronig 
relation of Purcell (1969; see eq.E]) if we assume that the extinction is caused by silicate dust 
and carbon dust 


H 


1 


A 


int 


.Csil 


(14) 


1.086 X 37r2 H 

where Fc and Fgn are respectively the dimensionless factor F seen in eq.[3]for carbon dust 
and silicate dust. For moderately elongated grains, Fsa ~ 0.7 for silicate materials and 
Fc ~ 1.25 for conducting, graphitic materials (see Figure [T]). For the selected 16 sources, 
the extinction per H nucleon Ax/N^i is only known over a limited range of wavelengths. 
Therefore, we hrst use eq.[9] to extrapolate the UV extinction of Valencic et ah (2004) to 
A = 912 A@ For 0.001 < A“^ < 0.3 pm”^, we adopt the model Ax/N^i values of WDOl 
and Li & Draine (2001b) for the diffuse ISM. This is justihed since the near- and mid- 
IR extinction at A > 0.9 pm does not seem to vary much among different environments 
(see Draine 2003a, Wang & Jiang 2014, Wang et al. 2013, 2014). In this way, we obtain 


R(nt = / Ax/N-^dX for each source (see Tabled]). Since < Rjnt, the volume of 
carbon dust per H nucleon (Vc/H) derived from eq.dHis a lower limit if we replace Hint by 
A[^^. In the following, unless otherwise stated, by Hint we really mean A[ 


int * 


What is the silicate dust volume per H nucleon 14ii/H? As shown in §11 this depends 
on the assumed interstellar Fe, Mg and Si abundances. In the literature, the proto-Sun 
abundances (Fodders 2003), the solar abundances (Asplund et al. 2009), the abundances of 
young F and G stars (Soha & Meyer 2001) and/or unevolved early B-type stars (Przybilla 
et al. 2008) have been adopted to represent the interstellar standard abundances. In this 
work we will only consider the proto-Sun abundances and the abundances of early B stars 
since the solar abundances and the abundances of young F and G stars are intermediate 
between that of the proto-Sun and B stars (see Footnote-[T|) . The proto-Sun abundances of 
Fe, Mg and Si result in 14ii/H 3.17 x 10“^^cm^H“^ while the B-star abundances lead to 
Rsii/H 2.52 X 10“^^cm^H“^. Assuming the carbon dust to be graphite we determine 
Vc/H from eq.dUfor each sightline. The required carbon depletion is 

[C/H]d„,t^(Hc/H)pc/12mH , (15) 


^Gordon et al. (2009) studied the extinction curves of 75 Galactic sightlines obtained with FUSE at 
905 < A < 1187 A and the International Ultraviolet Explorer (lUE) at 1150 < A < 3300 A. They found that 
the extrapolation of the UV extinction at 3.3 < A“^ < 8.7 obtained by lUE is generally consistent 

with the far-UV extinction at 8.4 < A“^ < 11 /im“^ obtained by FUSE. 

^°As shown in Figure [T] for moderately elongated or flattened spheroids the Fc factors are essentially the 
same for graphite and amorphous carbon. Therefore, the results derived from eq.[T4]also hold for amorphous 
carbon. 
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where pc is the mass density of the carbon dnst [pc ~ 2.24gem ^ for graphite). 

In Tabled] we list the total carbon volnme Vc/H and carbon depletion per H 

atom reqnired to acconnt for the observed extinction for each sightline. For a given set of 
interstellar Fe, Mg and Si reference abnndances, the nncertainties on mainly arise 

from Fc, Fsii and [Si/Hjjgj^. For moderately elongated or flattened spheroids, the Fq and 
Fsii valnes adopted here can vary at most ~20%. For a given reference standard (proto-Snn 
or B stars), the nncertainty on [Si/Hjjgj^ is within ~5%. Therefore, the nncertainties on 
are within ~25%. We note that, with Ay ~ 2.16 ± 0.04 mag (Gordon et ah 2009) 
and Nh ~ 1.78 x lO^^cm'^H (Cartledge et al. 2004), the line of sight toward HD 206773 is 
very dnsty in spite it has a large gas-phase carbon abnndance ([C/Hjgg^g 464 ± 57ppm; 
Parvathi et al. 2012). With Ay/N^i ~ 1.12 x 10“^^ mag cm^ the HD 206773 sightline 
is twice as dnsty as the Galactic diffnse ISM of which Ay/N^i ~ 5.3 x 10“^^ mag cm^ H~^. 
Therefore, one reqnires an nnnsnally large carbon depletion ([G/Hj^^^^ ^ 733ppm if the 
proto-Snn Fe, Mg, and Si abnndances are adopted, or [G/Hj^^g^ 774 ppm if the B-star Fe, 

Mg, and Si abnndances are adopted), mnch larger than that of the other sightlines considered 
in this work. In the following, nnless otherwise stated, we will not inclnde HD 206773 for 
correlation stndies. 

In Fignre H] we plot the carbon depletion [G/Hj^^^^ as a fnnetion of Cg, the strength 
of the 2175 A extinction bnmp. It can be seen from Fignre 0] that [G/Hj^^g^. shows an 
increasing trend with Cg, with a Pearson linear correlation coefficient of F ~ 0.60, and a 
Kendall correlation coefficient of r ~ 0.39 and a corresponding probability p ~ 0.04 of a 
chance correlation]^ This is independent of the assnmed interstellar Fe, Mg and Si reference 
abnndances which wonld only affect the intercept. In Fignre [5] we plot the carbon depletion 
[G/H]j^g^ as a fnnetion of the strength of the nonlinear far-UV extinction rise. Fignre [5] 
shows no detectable correlation between [G/Hj^^^^ and C 4 for which the Pearson correlation 
coefficient is i? ~ —0.20 and the Kendall r is ~ —0.01 and p ~ 0.96. Again, these correlation 
coefficients are independent of the assnmed interstellar Fe, Mg and Si reference abnndances. 


3. Carbon Depletion Inferred from the Extinction Modeling 

The carbon depletion [G/Hj^^^^ derived from the Kramers-Kronig relation is independent 
of any dnst models except we jnst need to assnme that the observed extinction is cansed by 


set the level of confidence to be 95% (i.e., Kendall’s a = 0.05). 
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silicate dust and carbon dusto 

We now would like to derive the carbon depletion from fitting the observed extinction 
curve for each sightline. We consider the silicate-graphite interstellar grain model which 
consists of two separate dust components: amorphous silicate and graphite (Mathis et ah 
1977, Draine & Lee 1984). We adopt an exponentially-cutoff power-law size distribution for 
both components: drii/da = nE_Bia~°‘^ exp {—a/ac,i) for the size range of 50 A < a < 2.5/im, 
where a is the spherical radius of the dust, riu is the number density of H nuclei, drii is 
the number density of dust of type i with radii in the interval [a, a + da] , o;* and Qc^i are 
respectively the power index and exponential cutoff size for dust of type z, and Bi is the 
constant related to the total amount of dust of type i. The total extinction per H column 
at wavelength A is given by 

/ I dfi 

da —(16) 
uh aa 

where the summation is over the two grain types (i.e., silicate and graphite), Nn = f rindl 
is the H column density which is the H number density integrated over the line of sight I, 
and Cext,i(a, A) is the extinction cross section of grain type i of size a at wavelength A which 
can be calculated from Mie theory (Bohren & Huffman 1983) using the dielectric functions 
of “astronomical” silicate and graphite of Draine & Lee (1984). 

In htting the extinction curve, for a given sightline, we have six parameters: the size 
distribution power indices as and ac for silicate and graphite, respectively; the exponential 
cutoff sizes ac,s and ac,c, respectively; and Bs and Bq- We derive the silicon and carbon 
depletions from 

[Si/H]j^g^ = (nHSs/172mH) j da (47r/3) exp (-a/ac,s) , (17) 

[C^/H]dust = (^H5c/12mH)yda(47r/3)a^Pgi.aa"“°exp(-a/ac,c) , (18) 

where we assume a stoichiometric composition of MgFeSiOd for amorphous silicate. 

For a given sightline, we seek the best £t to the extinction between 0.3 /zm“^ and 8 
by varying the size distribution power indices as and ac, and the upper cutoff size parameters 
ac,s and ac,c- Following WDOl, we evaluate the extinction at 100 wavelengths Aj, equally 


^^The correlation coefficients derived in ij2]and shown in Figure |4] remain the same for both graphite and 
amorphous carbon. The required carbon depletion would decrease by a factor of Pgra/pac, where 

Pac ~ 1.8gcm“^ is the mass density of amorphous carbon. 
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spaced in In A. We use the Levenberg-Marquardt method (Press et al. 1992) to minimize 
which gives the error in the extinction £t: 


= 




(In A, 


obs 


\nA 


mod ) 


(Jt 


(19) 


where Aobs(Aj) is the observed extinction at wavelength A,, Ainod(Ai) is the extinction com¬ 
puted for the model at wavelength A* (see eq.[T6]), and the Oi are weights. Following WDOl, 
we take the weights = 1 for 1.1 < A“^ < 8 /im”^ and = 1/3 for A“^ <1.1 /im”^. 

In Figures [2],[3] we show the model fits. It can be seen from these figures that a simple 
mixture of silicate and graphite closely reproduces the observed UV/optical/near-IR extinc¬ 
tion of all 16 sightlines. The model parameters are tabulated in Table [2l In Figure [6] we 
examine the correlations between the strength of the 2175 A extinction bump (cg) with the 
carbon ([C/Hj^^^J and silicon ([Si/II]^^^^) depletions derived from htting the observed ex¬ 
tinction (see eqs. fTTlfTSj) . With a Pearson correlation coefficient of i? ~ 0.83 and a Kendall 
r ^ 0.49 and p ~ 0.01, we see that the carbon depletion correlates with the 2175 A bump, 
while the silicon depletion does not correlate with the 2175 A bump. We have also explored 
the relations between the strength of the nonlinear far-UV extinction rise {c'^) and the carbon 
and silicon depletions. As shown in Figure [71 we do not see any correlation. 


4. Discussion 

Interstellar dust plays an important role in various physical and chemical processes 
that take place in the ISM. Dust extinction remains the most direct way to investigate the 
properties of interstellar dust since Trumpler (1930) hrst established the existence of solid 
particles in interstellar space through the discovery of color excess. However, practically most 
of the extinction modeling efforts have been so far directed to the Galactic average extinction 
curve which is obtained by averaging over many clouds of different gas and dust properties. 
An investigation of dust extinction toward individual clouds is of special importance. 

In this work we explore the extinction and dust depletion of 16 individual sightlines of 
which Ry ranges from ~2.4 to ~4.3 and the average hydrogen density (rin) ranges over 
nearly three orders of magnitudes (see Figure [8]). The hydrogen column density and the gas- 
phase carbon abundance are known for these sightlines, with the latter determined from the 
strong 1334 A CII absorption line (Parvathi et al. 2012, Soha et al. 2011). Unlike Cardelli 
et al. (1996) who derived the interstellar gas-phase carbon abundance to be [C/H]gg^g 
(140 ± 20) ppm and found that [C/Hjgg^g shows no dependence on the physical condition of 
the gas, for the sightlines studied by Parvathi et al. (2012) and Sofia et al. (2011), the gas- 
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phase [C/H]gj^g appears to show a decreasiM trend for increasing valnes of (nn) (see FignreEl 
also see Figure 3 in Parvathi et ah 2012)0 

The lack of any statistically signihcant variation in the observed gas-phase [0/11]^^^^ 
among the six sightlines studied by Cardelli et ah (1996) is in strong contrast to the appre¬ 
ciable variation of [C/H]g^g with (uh) seen in the sightlines studied in this work. As shown in 
FigurelH the solid-phase silicon ([Si/H]^^gJ and carbon ([C/H]j^g^) abundances derived from 
fitting the observed extinction also exhibit substantial variations with (uh). The solid-phase 
carbon abundance ([C/H]j^j,J determined from the wavelength-integrated extinction (Aint) 
based on the Kramers-Kronig relation varies by a factor of >10 among these 16 sightlines. 
Therefore, the total abundance of carbon, the sum of the components residing in the gas 
and in the dust, varies signihcantly among these sightlines. This raises a question: is there 
indeed a reference standard for the interstellar carbon abundance? 

As shown in Table [H the analysis of the wavelength-integrated extinction based on the 
Kramers-Kronig relation shows that five of the 16 sources require in total a carbon abundance 
(i.e., [C/H]j^^^ + [C/H]g^) exceeding 500 ppm, much higher than that of the young F and G 
stars ([C/H]^ ^ 355 ±82 ppm, Sofia & Meyer 2001), the proto-Sun ([C/Hj^ 288 ±27ppm. 

Fodders 2003), the Sun ([C/H]q 269 ±31 ppm, Asplund et al. 2009), and the early B stars 

([C/H]^ 214 ± 20 ppm, Przybilla et al. 2008), which are often taken to be the interstellar 

carbon reference abundance. This is also true for the carbon depletion derived from htting 
the observed extinction with the silicate-graphite model. Mathis (1996) argued that a dust 
model may require fewer C, Si, Mg and Fe atoms to account for the observed extinction if 
the dust has a fluffy, porous structure since fluffy grains are more effective in absorbing and 
scattering the UV/optical starlight than compact grains (on a per unit mass basis). However, 
as demonstrated in Li (2005), fluffy dust is not able to appreciably reduce the consumption 
of dust-forming atoms. As can be seen in eq.El although the dust volume Vdust does increase 
when the dust has voids and becomes porous, the F-factor decreases and therefore, for a given 
amount of dust material, the wavelength-integrated extinction J Ax/N^^dX does not change 
much (see Figures 1,2 in Li 2005). Alternatively, one may argue that the gas-phase [C/H]gg^g 
abundances derived from the strong CII absorption line at A ~ 1334 A may be overestimated 
since accurate measures of the absorption in the wings of this line is very difficult. However, 
the weak CII line at A ~ 2325 A gives even higher [C/H]gj^g abundances (Soha et al. 2004, 
Parvathi et al. 2010). Therefore, it seems real that the total C/H abundance of some of our 
16 sigtlines is well above that of the popular interstellar reference standards. 


^^Cardelli et al. (1996) determined [C/H]g^g from the weak intersystem line of CII] at 2325 A obtained 
with the HST Goddard High Resolution Spectrograph (GHRS) for six sightlines toward stars within 600 pc 
of the Sun. These sightlines include examples exhibiting a wide range of extinction variation. 
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In fitting the observed extinction (see ^ and Table [2]), seven out of the 16 sources 
require [Si/Hj^^^^ to exceed the proto-Sun silicon abundance by >50%, three of them are 
within ~ 20% of the proto-Sun silicon abundance, and six of them require less silicon than 
that of the proto-Sun. The variations in the silicon and carbon depletions toward differ¬ 
ent sightlines imply various physical processing (e.g., growth, coagulation, destruction, and 
photo-processing) of interstellar dust and variations in the composition within the Milky 
Way (Parvathi et ah 2012). 


Parvathi et al. (2012) found no correlation between the carbon depletion and 

the strength of the 2175 A extinction bump (cg), while in this work we hnd corre¬ 

lates with Cg (see ^and The major difference between our approaches and Parvathi et 
al. (2012) is that they assumed a reference standard of [C/Hjjgj^ = 464 ppm for all sightlines, 
while the approaches taken by us which make use of the Kramers-Kronig relation (see §2]) 
or model the observed extinction (see ^ result in considerable variation in the total C/H 
abundance. We note that for the Kramers-Kronig-relation-based approach we have assumed 
a reference standard for Si, Mg and Fe (see One may argue that the correlation between 
and Cg derived from this approach may not be real: with [Si/Hj^^g^., [Mg/Hj^^g^ 
and [Fe/Hjj^g^ hxed at a reference standard, one might expect the required C/H depletion 
to increase with the extinction. However, this is not true. With [Si/Hj^^g^, [Mg/Hj^^g^. and 
[Fe/H]j^g^ hxed, we do expect [C/Hj^^g^ to increase with Ha/Wi, but we do not necessarily 
expect [C/H]^^g^ to increase with Cg since the latter is a quantity normalized by Ay and 

therefore the common correlation of [C/Hj^^^g^ and Cg with Ay is cancelled. This can be 

tested by hxing [C/Hjjgjyj and determining the Si/H depletion from the Kramers-Kronig re¬ 
lation and then examining the correlation between Cg and [Si/Hj^^g^. Assuming the same 
total C/H abundance of [C/H]jgj^ for all 16 sightlines, the carbon dust volume per H nucleon 
for a given sightline can be derived from its gas-phase C/H abundance: 


^ = {|C/H],s„-|C/H|j„}xl2mH/pc ■ (20) 

We use the Kramers-Kronig relation to determine the silicate volume per H nucleon from 
the wavelength-integrated extinction: 


Vsii _ f 1 ) f Aint 

H VFsiJ V 1.086 X 37 r 2 ^hJ 

Then we derive the required Si/H depletion from 


( 21 ) 


[Si/H]dugt ~ (14il/H) Psil//^sil^H 5 


( 22 ) 


where /Xsii = 172 is the molecular weight for MgFeSi 04 . Following Parvathi et al. (2012), 
we adopt [C/H]jgj^ = 464 ppm. In Figure [9] we plot [Si/H]^^g^ against Cg and c'^ and see no 
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correlation. This justifies the correlation between [C/H]^^g^ and Cg derived in the Kramers- 
Kronig-relation-based approach presented in ^ which assumes fixed [Si/Hj^^^g^, 
and [Fe/H]^^g^ depletions. 

We note that although the carbon abundance [C/H]jgj^ is underabundant in accounting 
for the observed extinction, no matter what reference standard (proto-Sun, Sun, early B stars, 
or young F and G stars) is adopted, the oxygen abundance remains largely unaccounted for 
(see Whittet 2010). Recently, Poteet et al. (2015) argued that the missing oxygen may reside 
in large, micrometer-sized H 2 O ice grains which are nearly opaque to IR radiation. Possible 
evidence for the presence of /im-sized interstellar grains have been found by Wang et al. 
(2014), who modeled the mid-IR extinction toward a variety of interstellar environments 
using a population of /xm-sized grains. Furthermore, the recent analysis of three possibly 
interstellar grains returned by the Stardust spacecraft also provides tentative evidence for 
the presence of /xm-sized grains in the ISM (Sterken et al. 2014, Westphal et al. 2014). These 
/xm-sized grains, although gray in the UV/optical/near-IR, will contribute to the wavelength- 
integrated extinction Rjnt and therefore reduce the consumption of silicon and carbon (i.e., 
decreasing and [Si/H]^^^J. 

The correlation between the 2175 A bump strength and the carbon depletion 
inferred both from the Kramers-Kronig relation and from fitting the observed extinction 
with a specific dust model supports the hypothesis of carbonaceous grains (e.g., graphite or 
PAHs) as the possible carriers of the 2175 A bump0 The lack of correlation between the 
silicon depletion and the 2175 A bump argues against small silicates or (Mg, Si) oxides as its 
carrier (see Duley 1985). 

Finally, we examine the relations between the carbon and silicon depletions with Ry^. 
As shown in Figure fTOl the carbon depletion appears to moderately correlate with Ry^. Ry is 
roughly indicative of the dust size: the smaller Ry is, the steeper is the far-UV extinction and 
the richer in small dust is the dust size distribution. This suggests that the far-UV extinction 
is more likely produced by small carbon dust than by small silicate dust, probably arising 
from the blue wing of the a-a* electronic absorption profile of graphitic materials. We note 
that Ry provides a better characterization of the far-UV extinction than C 4 as the latter 
only accounts for the far-UV nonlinear extinction rise while the linear background 
surely also contributes to the far-UV extinction. In Figure [11] we plot the C/H and Si/H 
depletions against 03 . As expected, we do not see any correlation because, similar to c'^, 


Amorphous carbon (Bussoletti et al. 1987), graphitized hydrogenated amorphous carbon (HAG; Mennella 
et al. 1996), and nano-sized HAG (Schnaiter et al. 1998) have also been proposed as candidate materials. 
However, they are not always successful in reproducing the observed invariant peak position and variable 
strength and width of the 2175 A bump. 
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alone does not characterize the far-UV extinction. 


5. Conclusion 

We have studied the extinction and dust depletion in 16 Galactic sightlines in order to 
probe the role of carbon in the UV extinction, with special attention paid to the 2175 A 
extinction bump. These sightlines have their UV/optical/near-IR extinction, gas-phase car¬ 
bon abundances, and hydrogen column densities well determined, therefore allowing us to 
quantitatively explore the relations between the silicon and carbon depletions and the UV 
extinction. Our principal results are as follows: 

1. We have derived the carbon depletion from the Kramers-Kronig relation 

which relates the wavelength-integrated extinction to the total dust volume. We hnd 
that [C/H]j^g^ intends to correlate with the strength of the 2175 A extinction bump. 
This is less model-dependent since it does not require the knowledge of the detailed 
optical properties and size distribution of the dust. 

2. We have derived the carbon [C/H]j^g^ and silicon [Si/Hj^^g^. depletions from htting 
the observed UV/optical/near-IR extinction with a mixture of silicate dust and car¬ 
bonaceous dust. We hnd that [C/H]^^g^ correlates with the 2175 A extinction bump 
and [Si/H]j^g^. does not correlate with the 2175 A bump. This supports graphite or 
PAH molecules as the possible carrier of the 2175 A bump, while argues against the 
hypothesis of silicate or (Mg, Si) oxides as its carrier. 

3. A moderate correlation between [C/H]^^g^ and l/Ry is derived from the extinction 
modeling. This suggests that the far-UV extinction is more likely produced by small 
carbon dust than by small silicate dust. 


We thank B.W. Jiang, Q. Li, B. Mills, J.Y. Seok, and Y.X. Xie for helpful discussions. 
We are particularly grateful to S. Federman, S. Wang, and the anonymous referee for their 
suggestions which signihcantly improved the quality of this paper. We are supported in part 
by NSF AST-1109039, NNX13AE63G, and NSFG 11173019, 11273022 and the University of 
Missouri Research Board. 


REFERENCES 


Allamandola, L.J., Tielens, A.G.G.M., & Barker, J.R. 1985, ApJ, 290, L25 



16 


Asplund, M., Grevesse, N., Sauval, A.J., & Scott, P. 2009, ARA&A, 47, 481 

Bohren, C.F., & Huffman, D.R. 1983, Absorption and Scattering of Light by Small Particles 
(New York: Wiley) 

Bussoletti, E., Colangeli, L., Borghesi, A., & Orofino, V. 1987, A&AS, 70, 257 

Cardelli, J.A., Clayton, G.C., & Mathis, J.S. 1989, ApJ, 345, 245 (CCM) 

Cartledge, S.I.B., Lauroesch, J.T., Meyer, D.M., & Soha, U.J. 2004, ApJ, 613, 1037 

Cartledge, S.I.B., Lauroesch, J.T., Meyer, D.M., & Soha, U.J. 2006, ApJ, 641, 327 

Cecchi-Pestellini, C., Malloci, G., Mulas, G., Joblin, C., & Williams, D. A. 2008, A&A, 486, 
L25 

Chiar, J. E., Tielens, A. G. G. M., Adamson, A. J., & Ricca, A. 2013, ApJ, 770, 78 

Davis, A., Leigh, C.R.A., Hampton, R.N., Saunders, G.A., & Parker, S.C. 1988, J. Mater. 
Sci. Lett., 7, 415 

Draine, B.T., & Lee, H.M. 1984, ApJ, 258, 89D 
Draine, B.T., & Malhotra, S. 1993, ApJ, 414, 632 
Draine, B.T. 2003a, ARA&A, 41, 241 

Draine, B.T. 2003b, in The Cold Universe, Saas-Fee Advanced Course Vol. 32, ed. D. Pfen- 
niger (Berlin; Springer-Verlag), 213 

Duley, W.W., Jones, A.P., & Williams, D.A. 1989, MNRAS, 236, 709 

Duley, W.W. 1985, Ap&SS, 112, 321 

Fitzpatrick, E.L., & Massa, D. 1990, ApJS, 72, 163 

Furton, D. G., Laiho, J. W., & Witt, A. N. 1999, ApJ, 526, 752 

Gordon, K.D., Cartledge, S., & Clayton, G.C. 2009, ApJ, 705, 1320 

Greenberg, J.M., & Li, A. 1996, A&A, 309, 258 

Henning, Th. 2010, ARA&A, 48, 21 

Hildebrand, R.H., & Dragovan, M. 1995, ApJ, 450, 663 



17 


Jenkins, E.B. 2009, ApJ, 700, 1299 

Joblin, C., Leger, A., & Martin, P. 1992, ApJ, 393, L79 

Jones, A.P., Fanciullo, L., Kohler, M., et al. 2013, A&A, 558, AA62 

Kemper, F., Vriend, W.J., & Tielens, A.G.G.M. 2004, ApJ, 609, 826 

Kim, S.H., Martin, P.G., & Hendry, P.D. 1994, ApJ, 422, 164 

Lee, H.M., & Draine, B.T. 1985, ApJ, 290, 211 

L%er, A., & Puget, J. 1984, A&A, 137, L5 

Li, A., & Greenberg, J.M. 1997, A&A, 323, 566 

Li, A., & Draine, B.T. 2001a, ApJ, 550, L213 

Li, A., & Draine, B.T. 2001b, ApJ, 554, 778 

Li, A. 2005, ApJ, 622, 965 

Li, M.P., Shi, Q.J., & Li, A. 2008, MNRAS, 391, L49 
Lodders, K. 2003, ApJ, 591, 1220 

Mathis, J.S., Rumpl, W., & Nordsieck, K.H. 1977, ApJ, 217, 425 
Mathis, J.S. 1994, ApJ, 422, 176 
Mathis, J. S. 1996, ApJ, 472, 643 

Mennella, V., Golangeli, L., Palumbo, P., et al. 1996, ApJ, 464, L191 

Mulas, G., Zonca, A., Gasu, S., & Gecchi-Pestellini, G. 2013, ApJS, 207, 7 

Parvathi, V.S., Soha, U.J., Murthy, J., & Babu. B.R.S. 2012, ApJ, 760, 36 

Parvathi, V. S., Soha, U. J., Murthy, J., & Babu, B. R. S. 2010, Bull. Astron. Soc. India, 
38, 93 

Pendleton, Y. J., & Allamandola, L. J. 2002, ApJS, 138, 75 
Poteet, G. A., Whittet, D. G. B., & Draine, B. T. 2015, ApJ, 801, 110 



18 


Press, W.H., Teukolsky, S.A., Vetterling, W.T., & Flannery, B.P. 1992, Numerical Recipes 
in FORTRAN: The Art of Scientific Computing (2nd ed.; Cambridge: Cambridge 
Univ. Press) 

Przybilla, N., Nieva, M.-F., & Butler, K. 2008, ApJ, 688, L103 
Purcell, E.M. 1969, ApJ, 158, 433 

Schnaiter, M., Mutschke, H., Dorschner, J., Henning, T., & Salama, F. 1998, ApJ, 498, 486 
Snow, T.P., & Witt, A.N. 1996, ApJ, 468, L65 
Sofia, U.J., & Meyer, D.M. 2001, ApJ, 558, L147 

Sofia, U.J. 2004, in ASP Conf. Ser. 309, Astrophysics of Dust, ed. A.N. Witt, G.C. Clayton, 
& B.T. Draine (San Francisco: ASP), 393 

Sofia, U. J., Lauroesch, J. T., Meyer, D. M., & Cartledge, S. I. B. 2004, ApJ, 605, 272 
Soha, U.J., Parvathi, V.S., Babu, B.R.S., & Murthy, J. 2011, AJ, 141, 22 
Stecher, T.P., & Bonn, B. 1965, ApJ, 142, 1681 
Steglich, M., Jager, C., Rouille, G., et ah 2010, ApJ, 712, L16 

Sterken, V. J., Westphal, A. J., Altobelli, N., et al. 2014, Meteorit. Planet. Sci., 49, 1680 

Trumpler, R.J. 1930, PASP, 42, 214 

Wang, S., & Jiang, B.W. 2014, ApJ, 788, L12 

Wang, S., Li, A., & Jiang, B.W. 2014, Planet. Space Sci., 100, 32 

Wang, S., Gao, J., Jiang, B.W., Li, A., & Chen, Y. 2013, ApJ, 773, 30 

Weingartner, J.C., & Draine, B.T. 2001, ApJ, 548, 296 (WDOl) 

Westphal, A.J., et al. 2014, Science, 345, 786 

Whittet, D.C.B. 2003, Dust in the Galactic Environment (2nd ed; Bristol: loP) 

Whittet, D. C. B. 2010, ApJ, 710, 1009 

Witt, A. N.,& Vijh, U. P. 2004, in ASP Conf. Ser. 309, Astrophysics of Dust, ed. A.N. Witt, 
G.C. Clayton, & B.T. Draine (San Francisco: ASP), 115 



Valencic, L.A., Clayton, G.C., Sz Gordon, K.D. 2004, ApJ, 616, 912 
Zhang, K., Jiang, B. W., & Li, A. 2009, MNRAS, 396, 1247 



Fig. 1.— The F{eo; shape) factor as a function of the axial ratio a/b for silicate (solid line), Fe 203 (dotted 
line), amorphous carbon (dot-dashed line), and graphite (dashed line) of which the static dielectric constants 
are approximately eo~6-4, 16, 160 and oo (Li 2005). The grains are taken to be spheroidal with a and b 
being the semiaxis along and perpendicular to the symmetry axis of the spheroid, respectively. The blue 
lines are for oblates {a/b < 1) and the red lines are for prelates {a/b > 1). For modestly elongated {a/b < 3) 
or flattened {b/a < 3) silicate dust, the F factor is always smaller than unity. 


This preprint was prepared with the AAS lAERX macros v5.2. 
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Fig. 2.— Observed and model extinction curves of HD 027778, HD 037903, HD 157857, HD 185418, 
HD 192639, HD 198781, HD 206773, and HD 207198. The observed extinction curves are represented by 
the FM90 parametrization at > 3.3/im“^ and by the COM parametrization at < 3.3(open 
red diamonds). The filled black circles plot the optical data at D, B, V bands and the lUE data at 
3.3 < A“^ < 8.7 (taken from Gordon et al. 2009). The solid green line plots the model extinction 
curve which is a combination of silicate (dotted purple line) and graphite (dashed blue line). 
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Fig. 3. — Same as Figure[3]but for HD 122879, HD 147888, HD 091824, HD 116852, HD 203532, HD 210809, 
HD 232532, and HD 037061. For the latter six sources, the optical and lUE data are not shown. 
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Fig. 4. — Correlation diagrams between the carbon depletion with the strength of the 2175 A 

extinction bump (C 3 ). The carbon depletion is derived from the Kramers-Kronig relation (see eqs. Em, with 
the interstellar Fe, Mg and Si abundances taken to be that of proto-Sun (a) or B stars (b). The assumed 
interstellar Fe, Mg and Si reference abundances would not affect the correlation coefficients but the intercept. 
R is the Pearson correlation coefficient. Also labelled are the Kendall’s r coefficient and the significance 
level p. 




Fig. 5.— Correlation diagrams between the carbon depletion with the strength of the far- 

UV nonlinear extinction rise (C 4 ). The carbon depletion is derived from the Kramers-Kronig relation (see 
eqs. Em, with the interstellar Fe, Mg and Si abundances taken to be that of the proto-Sun (a) or B stars 
(b). 










23 




Fig. 6 .— Correlation diagrams between the strength of the 2175 A extinction bump ( 03 ) with the carbon 
depletion (a) and silicon depletion [Si/Hj^^^^ (b) derived from fitting the extinction of each sightline 

with a mixture of silicate dust and graphite dust (see lj3]). 




Fig. 7.— Correlation diagrams between the strength of the nonlinear far-UV extinction rise (C 4 ) with the 
carbon depletion [C/Hj^^g^. (a) and silicon depletion (b) derived from fitting the extinction of each 

sightline with a mixture of silicate dust and graphite dust (see SJ3]). 
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Fig. 8.— Silicon ([Si/H]^^j,^.) and carbon depletions as well as gas-phase carbon abundance 

([C/H]g^g) plotted as a function of the average hydrogen density (nn). The [Si/Hj^^^g^. and depletions 

are derived from fitting the extinction of each sightline with a mixture of silicate dust and graphite dust (see 
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Fig. 9. — Correlation diagrams between the Si/H depletion [Si/Hj^^^gj with the strength of the extinction 
bump (cg; a) and with the far-UV extinction rise {c'^, b). The Si/H depletion is derived from the Kramers- 
Kronig relation, with the total C/H abundance fixed at [C/H]jgj^ = 464 ppm. 
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Fig. 10. — Correlation diagrams between Ry^ and the carbon depletion (a) and silicon depletion 

[Si/H]j^g^ (b) derived from fitting the extinction of each sightline with a mixture of silicate dust and graphite 
dust (see SjS]). 



160 n 

1 - ^^^^^^^- 


^ R = 0.01 

140- 

HD 206773* „ 

T = - 0.01; p = 0.96 - 


silicate-graphite model ' 

120- 

- 

i 100- 

- 


- 

» 80- 

- 

, _ 

• 

I 

o 

CD 

H/!S 

. . • • • 

'—' 

: 

40 - 

• : 


• * ^ 


• • • 

20- 

• 


0 “b-'-1-'-1-'-1-'-1-'- 

0.0 0.1 0.2 0.3 0.4 0.5 

C2 


Fig. 11.— Correlation diagrams between the linear background (c^ with the carbon depletion 

(a) and silicon depletion [Si/HJ^^g^. (b) derived from fitting the extinction of each sightline with a mixture 

of silicate dust and graphite dust. 









Table 1. Extinction Parameters and the Carbon Depletion Required to Account for the Observed Extinction 


Star 

Nh 

(10^"^ cm“^ H) 

Av 

(mag) 

Rv 

<=1 

*^2 

^3 

C4 

XQ 

7 

^frit 

(magcm^ H~^) 

Proto-Sun^ 

Vc/H‘ 


B Stars^^ 

Vc/H‘ 


HD 027778 

2.29±1.20 

1.09±0.03 

2.79±0.38 

0.87±0.71 

0.32±0.01 

1.31±0.10 

0.24±0.01 

4.59±0.04 

1.21±0.03 

1.39E-25 

1.69E-27 

190 

2.05E-27 

231 

HD 037061*^ 

5.37±1.23 

2.40±0.21 

4.29A0.21 

1.54±0.00 

O.OOiO.OO 

0.31±0.00 

0.05±0.00 

4.57±0.00 

0.90±0.00 

1.19E-25 

1.19E-27 

134 

1.55E-27 

175 

HD 037903 

2.88±1.15 

1.49±0.03 

4.11A0.43 

1.17±0.06 

O.lliO.Ol 

1.04±0.09 

0.14±0.01 

4.60±0.04 

1.19±0.03 

1.46E-25 

1.86E-27 

209 

2.22E-27 

250 

HD 091824 

1.44±1.12 

0.60±0.06 

3.10±0.31 

0.74±0.41 

0.29±0.07 

1.19±0.35 

0.14±0.02 

4.62±0.09 

l.OOiO.lO 

1.16E-25 

l.llE-27 

125 

1.47E-27 

166 

HD 116852'^ 

1.05±1.20 

0.51±0.12 

2.42±0.37 

0.52±0.25 

0.38±0.10 

0.63±0.17 

O.OliO.Ol 

4.55±0.04 

0.78±0.07 

1.13E-25 

1.03E-27 

117 

1.40E-27 

158 

HD 122879 

2.19±1.26 

1.41±0.04 

3.17±0.20 

1.15±0.08 

0.15±0.02 

0.53±0.06 

0.09±0.01 

4.57±0.04 

0.74±0.02 

1.50E-25 

1.96E-27 

220 

2.32E-27 

261 

HD 147888 

5.89±1.20 

1.97±0.03 

4.08±0.18 

1.43±0.05 

0.03±0.01 

0.74±0.07 

0.07±0.01 

4.58±0.03 

0.94±0.02 

9.51E-26 

5.93E-28 

66.7 

9.54E-28 

108 

HD 157857 

2.75±1.17 

1.37±0.05 

3.15±0.30 

0.78±0.09 

0.22±0.02 

1.47±0.12 

O.lliO.Ol 

4.54±0.03 

0.96±0.02 

1.27E-25 

1.39E-27 

156 

1.75E-27 

197 

HD 185418 

2.57±1.17 

1.39±0.04 

2.67±0.20 

1.37±0.09 

0.16±0.02 

1.89±0.12 

0.14±0.01 

4.59±0.03 

1.03±0.02 

1.55E-25 

2.08E-27 

234 

2.45E-27 

275 

HD 192639 

3.02±1.17 

2.06±0.04 

3.12±0.16 

0.78±0.05 

0.27±0.01 

1.23±0.07 

0.06±0.01 

4.57±0.02 

0.96±0.02 

1.80E-25 

2.71E-27 

304 

3.07E-27 

345 

HD 198781 

1.41±1.15 

1.03±0.03 

2.77±0.29 

0.44±0.09 

0.42±0.02 

2.09±0.23 

0.13±0.02 

4.69±0.06 

1.36±0.03 

2.13E-25 

3.53E-27 

397 

3.89E-27 

438 

HD 203532*^ 

2.75±1.17 

0.94±0.11 

3.37±0.24 

0.76±0.16 

0.27±0.03 

1.50±0.25 

0.20±0.04 

4.59±0.01 

1.26±0.04 

1.57E-25 

2.13E-27 

240 

2.50E-27 

281 

HD 206773 

1.78±1.12 

2.16±0.04 

4.08±0.29 

0.84±0.07 

O.lSiO.Ol 

0.75±0.07 

0.04±0.01 

4.57±0.04 

1.04±0.02 

3.34E-25 

6.52E-27 

733 

6.88E-27 

774 

HD 207198 

4.79±1.23 

1.54±0.03 

2.68±0.11 

0.73±0.05 

0.35±0.01 

1.16±0.07 

0.20±0.01 

4.62±0.03 

1.04±0.02 

9.30E-26 

5.41E-28 

60.8 

9.02E-28 

102 

HD 210809 

2.14±1.20 

0.98±0.05 

3.10±0.31 

0.60±0.31 

0.32±0.02 

0.83±0.08 

0.15±0.01 

4.60±0.09 

1.GOTO.10 

8.24E-26 

2.77E-28 

31.1 

6.38E-28 

72.1 

HD 232522*^ 

1.55±1.12 

0.82±0.16 

3.05±0.41 

0.59±0.12 

0.38±0.07 

1.06±0.22 

0.23±0.06 

4.55±0.01 

0.93±0.03 

1.54E-25 

2.06E-27 

232 

2.42E-27 

273 


The wavelength-integrated extinction Aj^t — 


^ The interstellar Si, Mg and Fe abundances are taken to be that of proto-Sun (Lodders 2003) for which the total silicate volume per H atom is V^ii/H PS 3.17 X 10”^'^ cm^ 
The interstellar Si, Mg and Fe abundances are taken to be that of unevolved early B stars (Przybilla et al. 2008) for which V^u/H ci 2.52 X 10”^"^ cm^ 

Data taken from Valencic et al. (2004). The others are taken from Gordon et al. (2009). 

® The total volume of carbon dust per H atom (in unit of cm^ H~^) required to account for the observed extinction. 

^ The carbon depletion [C/H]j^g^ (in unit of ppm) required to account for the observed extinction. 
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Table 2. Model Parameters for Fitting the UV/Optical/Near-IR Extinction with A 


Mixture of Silicate/Graphite Dust 


Star 

Av/Nn 

(10-22 jnagcm^ H~^) 

as 

<^C,S 

(fim) 

QC 

ac,C 

(fim) 


[C/H]g,, 

(ppm) 

[C/H]d,,t 

(ppm) 

[Si/H]du.t 

(ppm) 

HD 027778 

4.76±2.50 

-3.1 

0.14 

-3.0 

0.10 

0.83 

93±32 

207 

48.0 

HD 037061 

4.47±1.10 

-2.0 

0.17 

-2.4 

0.05 

1.05 

99±36 

80.7 

57.5 

HD 037903 

5.17±2.07 

-2.7 

0.23 

-3.4 

0.20 

0.47 

382±56 

168 

58.6 

HD 091824 

4.17±3.27 

-3.1 

0.24 

-3.8 

0.25 

3.19 

202±24 

141 

49.0 

HD 116852 

4.86±5.67 

-2.7 

0.08 

-3.4 

0.15 

1.54 

116±102 

177 

43.0 

HD 122879 

6.44±3.71 

-2.7 

0.18 

-3.4 

0.10 

6.19 

365±94 

184 

64.0 

HD 147888 

3.34±0.68 

-2.1 

0.16 

-3.5 

0.20 

2.28 

98±23 

101 

35.0 

HD 157857 

4.98±2.13 

-2.5 

0.08 

-3.9 

0.40 

0.14 

112±19 

225 

31.0 

HD 185418 

5.05±2.16 

-2.0 

0.06 

-3.7 

0.30 

4.14 

173±34 

310 

36.0 

HD 192639 

6.82±2.65 

-2.6 

0.10 

-3.5 

0.10 

3.16 

138±24 

227 

60.8 

HD 198781 

7.30±5.96 

-3.0 

0.10 

-3.1 

0.10 

0.12 

193±31 

313 

72.9 

HD 203532 

3.42±1.51 

-3.3 

0.21 

-3.2 

0.15 

0.87 

85±28 

165 

28.8 

HD 206773 

12.1±7.64 

-2.7 

0.20 

-3.1 

0.05 

3.40 

464±57 

206 

143 

HD 207198 

3.22±0.83 

-2.8 

0.06 

-2.6 

0.05 

3.49 

69±21 

168 

21.0 

HD 210809 

4.58±2.58 

-3.2 

0.19 

-3.2 

0.15 

4.12 

321±70 

118 

28.0 

HD 232522 

5.29±3.96 

-3.4 

0.29 

-3.0 

0.05 

4.12 

317±46 

183 

63.9 




yc^ 



